Experimental and computational investigation of the effect of Hsc70 structural variants on inhibiting amylin aggregation.

The misfolding and aggregation of human islet amyloid polypeptide (hIAPP), also known as amylin, have been implicated in the pathogenesis of type 2 diabetes (T2D). Heat shock proteins, specifically, heat shock cognate 70 (Hsc70), are molecular chaperones that protect against hIAPP misfolding and inhibits its aggregation. Nevertheless, there is an incomplete understanding of the mechanistic interactions between Hsc70 domains and hIAPP, thus limiting their potential therapeutic role in diabetes. This study investigates the inhibitory capacities of different Hsc70 variants, aiming to identify the structural determinants that strike a balance between efficacy and cytotoxicity. Our experimental findings demonstrate that the ATPase activity of Hsc70 is not a pivotal factor for inhibiting hIAPP misfolding. We underscore the significance of the C-terminal substrate-binding domain of Hsc70 in inhibiting hIAPP aggregation, emphasizing that the removal of the lid subdomain diminishes the inhibitory effect of Hsc70. Additionally, we employed atomistic discrete molecular dynamics simulations to gain deeper insights into the interaction between Hsc70 variants and hIAPP. Integrating both experimental and computational findings, we propose a mechanism by which Hsc70's interaction with hIAPP monomers disrupts protein-protein connections, primarily by shielding the ß-sheet edges of the Hsc70-ß-sandwich. The distinctive conformational dynamics of the alpha helices of Hsc70 potentially enhance hIAPP binding by obstructing the exposed edges of the ß-sandwich, particularly at the ß5-ß8 region along the alpha helix interface. This, in turn, inhibits fibril growth, and similar results were observed following hIAPP dimerization. Overall, this study elucidates the structural intricacies of Hsc70 crucial for impeding hIAPP aggregation, improving our understanding of the potential anti-aggregative properties of molecular chaperones in diabetes treatment.

Zinc-Epigallocatechin-3-gallate Network-Coated Nanocomposites against the Pathogenesis of Amyloid-Beta.

The aggregation of amyloid beta (Aß) is a hallmark of Alzheimer's disease (AD), a major cause of dementia and an unmet challenge in modern medicine. In this study, we constructed a biocompatible metal-phenolic network (MPN) comprising a polyphenol epigallocatechin gallate (EGCG) scaffold coordinated by physiological Zn(II). Upon adsorption onto gold nanoparticles, the MPN@AuNP nanoconstruct elicited a remarkable potency against the amyloid aggregation and toxicity of Aß in vitro. The superior performance of MPN@AuNP over EGCG@AuNP was attributed to the porosity and hence larger surface area of the MPN in comparison with that of EGCG alone. The atomic detail of Zn(II)-EGCG coordination was unraveled by density functional theory calculations and the structure and dynamics of Aß aggregation modulated by the MPN were further examined by discrete molecular dynamics simulations. As MPN@AuNP also displayed a robust capacity to cross a blood-brain barrier model through the paracellular pathway, and given the EGCG's function as an anti-amyloidosis and antioxidation agent, this MPN-based strategy may find application in regulating the broad AD pathology beyond protein aggregation inhibition.

Fuzzy supertertiary interactions within PSD-95 enable ligand binding.

The scaffold protein PSD-95 links postsynaptic receptors to sites of presynaptic neurotransmitter release. Flexible linkers between folded domains in PSD-95 enable a dynamic supertertiary structure. Interdomain interactions within the PSG supramodule, formed by PDZ3, SH3, and Guanylate Kinase domains, regulate PSD-95 activity. Here we combined discrete molecular dynamics and single molecule Förster resonance energy transfer (FRET) to characterize the PSG supramodule, with time resolution spanning picoseconds to seconds. We used a FRET network to measure distances in full-length PSD-95 and model the conformational ensemble. We found that PDZ3 samples two conformational basins, which we confirmed with disulfide mapping. To understand effects on activity, we measured binding of the synaptic adhesion protein neuroligin. We found that PSD-95 bound neuroligin well at physiological pH while truncated PDZ3 bound poorly. Our hybrid structural models reveal how the supertertiary context of PDZ3 enables recognition of this critical synaptic ligand.

Different Forms of Disorder in NMDA-Sensitive Glutamate Receptor Cytoplasmic Domains Are Associated with Differences in Condensate Formation.

The N-methyl-D-aspartate (NMDA)-sensitive glutamate receptor (NMDAR) helps assemble downstream signaling pathways through protein interactions within the postsynaptic density (PSD), which are mediated by its intracellular C-terminal domain (CTD). The most abundant NMDAR subunits in the brain are GluN2A and GluN2B, which are associated with a developmental switch in NMDAR composition. Previously, we used single molecule fluorescence resonance energy transfer (smFRET) to show that the GluN2B CTD contained an intrinsically disordered region with slow, hop-like conformational dynamics. The CTD from GluN2B also undergoes liquid-liquid phase separation (LLPS) with synaptic proteins. Here, we extend these observations to the GluN2A CTD. Sequence analysis showed that both subunits contain a form of intrinsic disorder classified as weak polyampholytes. However, only GluN2B contained matched patterning of arginine and aromatic residues, which are linked to LLPS. To examine the conformational distribution, we used discrete molecular dynamics (DMD), which revealed that GluN2A favors extended disordered states containing secondary structures while GluN2B favors disordered globular states. In contrast to GluN2B, smFRET measurements found that GluN2A lacked slow conformational dynamics. Thus, simulation and experiments found differences in the form of disorder. To understand how this affects protein interactions, we compared the ability of these two NMDAR isoforms to undergo LLPS. We found that GluN2B readily formed condensates with PSD-95 and SynGAP, while GluN2A failed to support LLPS and instead showed a propensity for colloidal aggregation. That GluN2A fails to support this same condensate formation suggests a developmental switch in LLPS propensity.

Integrative structural dynamics probing of the conformational heterogeneity in synaptosomal-associated protein 25.

SNAP-25 (synaptosomal-associated protein of 25 kDa) is a prototypical intrinsically disordered protein (IDP) that is unstructured by itself but forms coiled-coil helices in the SNARE complex. With high conformational heterogeneity, detailed structural dynamics of unbound SNAP-25 remain elusive. Here, we report an integrative method to probe the structural dynamics of SNAP-25 by combining replica-exchange discrete molecular dynamics (rxDMD) simulations and label-based experiments at ensemble and single-molecule levels. The rxDMD simulations systematically characterize the coil-to-molten globular transition and reconstruct structural ensemble consistent with prior ensemble experiments. Label-based experiments using Förster resonance energy transfer and double electron-electron resonance further probe the conformational dynamics of SNAP-25. Agreements between simulations and experiments under both ensemble and single-molecule conditions allow us to assign specific helix-coil transitions in SNAP-25 that occur in submillisecond timescales and potentially play a vital role in forming the SNARE complex. We expect that this integrative approach may help further our understanding of IDPs.

Probing the adsorption behavior and free energy landscape of single-stranded DNA oligonucleotides on single-layer MoS2with molecular dynamics.

Interfacing single-stranded DNA (ssDNA) with 2D transition metal dichalcogenides are important for numerous technological advancements. However, the molecular mechanism of this process, including the nature of intermolecular association and conformational details of the self-assembled hybrids is still not well understood. Here, atomistic molecular dynamics simulation is employed to study the distinct adsorption behavior of ssDNA on a single-layer MoS2in aqueous environment. The ssDNA sequences [T10, G10, A10, C10, U10, (GT)5, and (AC)5] are chosen on the basis that short ssDNA segments can undergo a spontaneous conformational change upon adsorption and allow efficient sampling of the conformational landscape. Differences in hybridization is attributed to the inherent molecular recognition ability of the bases. While the binding appears to be primarily driven by energetically favorable van der Waalsπ-stacking interactions, equilibrium structures are modulated by the ssDNA conformational changes. The poly-purines demonstrate two concurrently competingπ-stacking interactions: nucleobase-nucleobase (intramolecular) and nucleobase-MoS2(intermolecular). The poly-pyrimidines, on the other hand, reveal enhancedπ-stacking interactions, thereby maximizing the number of contacts. The results provide new molecular-level understanding of ssDNA adsorption on the MoS2surface and facilitate future studies in design of functional DNA/MoS2structure-based platforms for DNA sequencing, biosensing (optical, electrochemical, and electronic), and drug delivery.

Out-of-Equilibrium Biophysical Chemistry: The Case for Multidimensional, Integrated Single-Molecule Approaches.

Out-of-equilibrium processes are ubiquitous across living organisms and all structural hierarchies of life. At the molecular scale, out-of-equilibrium processes (for example, enzyme catalysis, gene regulation, and motor protein functions) cause biological macromolecules to sample an ensemble of conformations over a wide range of time scales. Quantifying and conceptualizing the structure-dynamics to function relationship is challenging because continuously evolving multidimensional energy landscapes are necessary to describe nonequilibrium biological processes in biological macromolecules. In this perspective, we explore the challenges associated with state-of-the-art experimental techniques to understanding biological macromolecular function. We argue that it is time to revisit how we probe and model functional out-of-equilibrium biomolecular dynamics. We suggest that developing integrated single-molecule multiparametric force-fluorescence instruments and using advanced molecular dynamics simulations to study out-of-equilibrium biomolecules will provide a path towards understanding the principles of and mechanisms behind the structure-dynamics to function paradigm in biological macromolecules.

Ensemble Switching Unveils a Kinetic Rheostat Mechanism of the Eukaryotic Thiamine Pyrophosphate Riboswitch.

Thiamine pyrophosphate (TPP) riboswitches regulate thiamine metabolism by inhibiting the translation of enzymes essential to thiamine synthesis pathways upon binding to thiamine pyrophosphate in cells across all domains of life. Recent work on the Arabidopsis thaliana TPP riboswitch suggests a multi-step TPP binding process involving multiple riboswitch configurational ensembles and that Mg2+ dependence underlies the mechanism of TPP recognition and subsequent transition to the expression-inhibiting state of the aptamer domain followed by changes in the expression platform. However, details of the relationship between TPP riboswitch conformational changes and interactions with TPP and Mg2+ ¬¬in the aptamer domain constituting this mechanism are unknown. Therefore, we integrated single-molecule multiparameter fluorescence and force spectroscopy with atomistic molecular dynamics simulations and found that conformational transitions within the aptamer domain's sensor helices associated with TPP and Mg2+ ligand binding occurred between at least five different ensembles on timescales ranging from µs to ms. These dynamics are orders of magnitude faster than the 10 second-timescale folding kinetics associated with expression-state switching in the switch sequence. Together, our results show that a TPP and Mg2+ dependent mechanism determines dynamic configurational state ensemble switching of the aptamer domain's sensor helices that regulates the stability of the switch helix, which ultimately may lead to the expression-inhibiting state of the riboswitch. Additionally, we propose that two pathways exist for ligand recognition and that this mechanism underlies a kinetic rheostat-like behavior of the Arabidopsis thaliana TPP riboswitch.

Molecular insights on the dynamic stability of peptide nucleic acid functionalized carbon and boron nitride nanotubes.

The strategic approaches to the design of self-assembled hybrids of biomolecular systems at the nanoscale such as deoxyribonucleic acid (DNA) with single-wall carbon nanotubes (CNTs) and their structural analog, boron nitride nanotubes (BNNTs), rely on understanding how biomolecules recognize and mediate intermolecular interactions with the nanomaterial's surface. In this paper, we consider peptide nucleic acid (PNA), which is a synthetic analog of DNA, and investigate its interaction with a zigzag CNT and BNNT of similar diameter. The results based on the molecular dynamics method find that PNA provides definitive contrasts in the adsorption on the tubular surface in aqueous solution: it prefers to wrap along the circumferential direction on a (11,0) CNT, whereas it binds along the axial direction adopting an extended configuration on a (11,0) BNNT. Moreover, gas-phase Monte Carlo simulations show a dependence of the nanotube diameter on the calculated adsorption energy, with BNNTs exhibiting higher adsorption energy compared to CNTs, and the largest-diameter (25,0) tubular configuration facilitates encapsulation of PNA rather than PNA being adsorbed on its sidewall. The results are expected to be of relevance in the design of novel PNA-based archetypal hybrid materials for nanoscale applications in health-related areas including biosensing.

Acetylation of Aß42 at Lysine 16 Disrupts Amyloid Formation.

The residue lysine 28 (K28) is known to form an important salt bridge that stabilizes the Aß amyloid structure, and acetylation of lysine 28 (K28Ac) slows the Aß42 fibrillization rate but does not affect fibril morphology. On the other hand, acetylation of lysine 16 (K16Ac) residue greatly diminishes the fibrillization property of Aß42 peptide and also affects its toxicity. This is due to the fact that lysine 16 acetylated amyloid beta peptide forms amorphous aggregates instead of amyloid fibrils. This is likely a result of increased hydrophobicity of the K16-A21 region due to K16 acetylation, as confirmed by molecular dynamic simulation studies. The calculated results show that the hydrophobic patches of aggregates from acetylated peptides were different when compared to wild-type (WT) peptide. K16Ac and double acetylated (KKAc) peptide aggregates show significantly higher cytotoxicity compared to the WT or K28Ac peptide aggregates alone. However, the heterogeneous mixture of WT and acetylated Aß42 peptide aggregates exhibited higher free radical formation as well as cytotoxicity, suggesting dynamic interactions between different species could be a critical contributor to Aß pathology.

Dynamics of carbohydrate strands in water and interactions with clay minerals: influence of pH, surface chemistry, and electrolytes.

Carbohydrate hydrogels are extensively used in pharmaceuticals and engineered biomaterials. Molecular conformations, assembly, and interactions of the carbohydrate strands with stabilizers such as clay minerals in aqueous solution are difficult to quantify in experiments and the hydrogel properties remain largely a result of trial-and-error studies. We analyzed the assembly of gellan gum in aqueous solution and interactions with dispersed clay minerals in all-atomic detail using molecular dynamics simulation, atomic force microscopy (AFM), and comparisons to earlier measurements. Gellan strands associate at low pH values of 2 and gradually disassemble to double strands with weak association of -0.4 kcal per mole carbohydrate ring as the pH values increases to 9. Ionization of the carbonic acid side groups in the backbone extends the chains and accelerates the conformational dynamics via rapidly changing intramolecular ion bridges. Gellan interactions with clay minerals depend on the strength of electric triple layers between clay, cations, and anionic polymer strands, as well as weaker hydrogen bonds along the edges, which are tunable as a function of the clay surface chemistry, local ionic strength, and pH values. Interaction energies range from -4 to +6 kcal per mol carbohydrate ring and were most favorable for electric triple layers with high charge mobility, which can be achieved for intermediate cation exchange capacity of the clay mineral and high pH values to increase ionization of the clay edges and of the polymer. The findings provide understanding and help control the dynamics and stabilization of carbohydrate hydrogels by clay minerals.

Dynamics of self-assembled cytosine nucleobases on graphene.

Molecular self-assembly of cytosine (C n ) bases on graphene was investigated using molecular dynamics methods. For free-standing C n bases, simulation conditions (gas versus aqueous) determine the nature of self-assembly; the bases prefer to aggregate in the gas phase and are stabilized by intermolecular H-bonds, while in the aqueous phase, the water molecules disrupt base-base interactions, which facilitate the formation of π-stacked domains. The substrate-induced effects, on the other hand, find the polarity and donor-acceptor sites of the bases to govern the assembly process. For example, in the gas phase, the assembly of C n bases on graphene displays short-range ordered linear arrays stabilized by the intermolecular H-bonds. In the aqueous phase, however, there are two distinct configurations for the C n bases assembly on graphene. For the first case corresponding to low surface coverage, the bases are dispersed on graphene and are isolated. The second configuration archetype is disordered linear arrays assembled with medium and high surface coverage. The simulation results establish the role of H-bonding, vdW π-stacking, and the influence of graphene surface towards the self-assembly. The ability to regulate the assembly into well-defined patterns can aid in the design of self-assembled nanostructures for the next-generation DNA based biosensors and nanoelectronic devices.

Theoretical study of gas and solvent phase stability and molecular adsorption of noncanonical guanine bases on graphene.

The gas and solvent phase stability of noncanonical (Gua)n nucleobases is investigated in the framework of dispersion-corrected density functional theory (DFT). The calculated results strongly support the high tendency for the dimerization of (Gua)n bases in both gas and solvent phases. An interplay between intermolecular and bifurcated H-bonds is suggested to govern the stability of (Gua)n bases which bears a correlation with the description of dispersion correction terms employed in the DFT calculations. For example, a higher polarity is predicted for (Gua)n bases by the dispersion-corrected DFT in contrast to the non-polar nature of (Gua)3 and (Gua)4 predicted by the hybrid meta-GGA calculations. This distinct variation becomes significant under physiological conditions as polar (Gua)n is likely to exhibit greater stabilization in the gas phase compared to solvated (Gua)n. Graphene acting as a substrate induces modification in base configurations via maximization of π-orbital overlap between the base and substrate. In solvent, the substrate-induced effects are further heightened with lowering of the dipole moments of (Gua)n as also displayed by the corresponding isosurface of the electrostatic potential. The graphene-induced stability in both gas and solvent phases appears to fulfill one of the prerequisite criteria for molecular self-assembly. The DFT results therefore provide atomistic insights into the stability and molecular assembly of free-standing noncanonical (Gua)n nucleobases which can be extended to understanding the self-assembly process of functional biomolecules on 2D materials for potential biosensing applications.

Hierarchical Self-Assembly of Noncanonical Guanine Nucleobases on Graphene.

Self-assembly characterizes the fundamental basis toward realizing the formation of highly ordered hierarchical heterostructures. A systematic approach toward the supramolecular self-assembly of free-standing guanine nucleobases and the role of graphene as a substrate in directing the monolayer assembly are investigated using the molecular dynamics simulation. We find that the free-standing bases in gas phase aggregate into clusters dominated by intermolecular H-bonds, whereas in solvent, substantial screening of intermolecular interactions results in π-stacked configurations. Interestingly, graphene facilitates the monolayer assembly of the bases mediated through the base-substrate π-π stacking. The bases assemble in a highly compact network in gas phase, whereas in solvent, a high degree of immobilization is attributed to the disruption of intermolecular interactions. Graphene-induced stabilization/aggregation of free-standing guanine bases appears as one of the prerequisites governing molecular ordering and assembly at the solid/liquid interface. The results demonstrate an interplay between intermolecular and π-stacking interactions, central to the molecular recognition, aggregation dynamics, and patterned growth of functional molecules on two-dimensional nanomaterials.

Ab initio study on the noncovalent adsorption of camptothecin anticancer drug onto graphene, defect modified graphene and graphene oxide.

The application of graphene and related nanomaterials like boron nitride (BN) nanosheets, BN-graphene hybrid nanomaterials, and graphene oxide (GO) for adsorption of anticancer chemotherapeutic camptothecin (CPT) along with the effect on electronic properties prior to functionalization and after functionalization has been reported using density functional theory (DFT) calculations. The inclusion of dispersion correction to DFT is instrumental in accounting for van der Waals π-π stacking between CPT and the nanomaterial. The adsorption of CPT exhibits significant strain within the nanosheets and noncovalent adsorption of CPT is thermodynamically favoured onto the nanosheets. In case of GO, surface incorporation of functional groups result in significant crumpling along the basal plane and the interaction is basically mediated by H-bonding rather than π-π stacking. Docking studies predict the plausible binding of CPT, CPT functionalized graphene and GO with topoisomerase I (top 1) signifying that CPT interacts through π stacking with AT and GC base pairs of DNA and in presence of nano support, DNA bases preferentially gets bound to the basal plane of graphene and GO rather than the edges. At a theoretical level of understanding, our studies point out the noncovalent interaction of CPT with graphene based nanomaterials and GO for loading and delivery of anticancer chemotherapeutic along with active binding to Top1 protein.

Adsorption of 3,4-dihydroxybenzoic acid onto hematite surface in aqueous medium: importance of position of phenolic -OH groups and understanding of the same using catechol as an auxiliary model.

A comparative adsorption kinetics, isotherms, dissolution and surface complexation of 3,4-dihydroxybenzoic acid (3,4-DHBA) and 1,2-dihydroxybenzene (catechol) at the hematite/electrolyte interface were investigated. The kinetics at pH 10 and 298.15K suggested that the adsorption behaviour of 3,4-DHBA and catechol onto hematite surface is similar and attain same equilibration time of 60 min. The adsorption kinetics data of 3,4-DHBA and catechol fit the pseudo-second-order kinetic equation of nonlinear form best. The adsorption density of 3,4-DHBA at pH≥9 increases and thereby mimics the behaviour of catechol. The solubility of hematite depends on both pH of the suspension and concentration of adsorbate. The inner-sphere complex is formed by 3,4-DHBA and catechol onto hematite surface but the mode orientation is likely to be different in the pH range 5-8 and 9-10. The advance microscopic scanning in conjunction with the vibration spectroscopy would provide better pictorial presentation of the mode of orientation of 3,4-DHBA and catechol onto hematite surface at different pH.

Density functional and molecular docking studies towards investigating the role of single-wall carbon nanotubes as nanocarrier for loading and delivery of pyrazinamide antitubercular drug onto pncA protein.

The potential biomedical application of carbon nanotubes (CNTs) pertinent to drug delivery is highly manifested considering the remarkable electronic and structural properties exhibited by CNT. To simulate the interaction of nanomaterials with biomolecular systems, we have performed density functional calculations on the interaction of pyrazinamide (PZA) drug with functionalized single-wall CNT (fSWCNT) as a function of nanotube chirality and length using two different approaches of covalent functionalization, followed by docking simulation of fSWCNT with pncA protein. The functionalization of pristine SWCNT facilitates in enhancing the reactivity of the nanotubes and formation of such type of nanotube-drug conjugate is thermodynamically feasible. Docking studies predict the plausible binding mechanism and suggests that PZA loaded fSWCNT facilitates in the target specific binding of PZA within the protein following a lock and key mechanism. Interestingly, no major structural deformation in the protein was observed after binding with CNT and the interaction between ligand and receptor is mainly hydrophobic in nature. We anticipate that these findings may provide new routes towards the drug delivery mechanism by CNTs with long term practical implications in tuberculosis chemotherapy.

Density functional study on the adsorption of the drug isoniazid onto pristine and B-doped single wall carbon nanotubes.

The current study explores a new strategy to incorporate single wall carbon nanotubes (SWNTs)/doped SWNTs as carrier modules in target-specific administration of antitubercular chemotherapeutics through covalent and noncovalent functionalization onto the nanotube sidewall. Density functional studies illustrate that noncovalent functionalization of isoniazid (INH) is preferred over covalent attachment, exhibiting low adsorption energy values, HOMO-LUMO gap and comparison of quantum molecular descriptors performed in (5,5) and (9,0) SWNT systems. Substitution doping of boron facilitates the adsorption of INH onto the otherwise inert nanotube. Frontier orbital analysis reveals reorientation of electronic charge in the nanotubes after functionalization, the effect being more pronounced in the case of doped nanotubes. The charge transfer is significant in covalent functionalization of INH via the B-dopant atom, whereas in noncovalent functionalization a small amount of charge transfer is noted. Solvation studies demonstrate the dissolution of INH in B-doped (5,5) and (9,0) SWNTs to be higher compared to pristine nanotube-INH complexes. Functionalization of nanotubes via covalent and noncovalent means can foster pioneering prospects especially for experimental studies in this area of research.